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Abstract 

All one-massless-loop Feynman diagrams could be written like a linear combination of scalar 
boxes, triangles an bubbles in n dimensions plus rational terms. However, the four-point scalar 
integrals in n + 2 dimensions are free of infrared divergences. We are going to change the dimensions 
of the scalar boxes n — > n + 2 and the using of this degree of freedom to simplify the computation 
of coefficients of the decomposition. 
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1 Introduction 



Since many years, one have try to calculate analytically Feynman Diagrams. The aim is the com- 
putation of amplitudes and cross sections. In fact, at the LHC, we want to discover new physics and 
new particles like Higgs by the interaction of two protons. The background is constituted by many 
QCD reactions. So the knowlegde of this pone is necessary, if we want to detect a new particle. But 
the amplitudes depend to an unphysical energy, and this dependance decreases with the order of the 
development. So to have a good prediction, we have to calculate each reaction at NLO. At this time, 
all 2 — > 2 and 2^3 processes are known at NLO, but it remains the 2 — » N, with N > 3. So we have 
to calculate one-loop diagrams with many ingoing legs. 

Forty years ago, Passarino and Veltman gave a first method to reduce diagrams. This method is 
not really efficient and doesn't use the mathematical symmetry of a loop, for example, unitarity. They 
showed that we can write a diagram like a linear combination of scalar integrals. Then since ten years, 
Bern and al. [U H] reduce diagrams thanks to unitarity. It simplifies computation. In 2004, Britto 
and al. El [7l [9] gave a very efficient method to calculate the coefficients in front of the four- 
point scalar integrals. Then, Mastrolia [10] found a way to express the coefficients in front of triangles. 
Recently Forde [TT] in one hand and Papadopoulos and al. [T2] in a other hand gave some algorithms 
to obtain directly the coefficients in front of the scalar massless boxes, triangles and bubbles. Finally 
Kilgore [T3] improves the Forde algorithms to a massive loop. But in this algorithm, less the scalar 
integrals have legs, more the coefficients are difficult to calculate because some free parameters appear. 
In fact, as they use unitarity and cuts, we need four cuts to define a loop momenta. But we have only 
three cuts (respectively two cuts) in a triangle (resp. in a bubble). So some free parameters remains. 

In this paper, I would like to give a way to eliminate this degree of freedom in front of the triangles 
in a massless loop. Just after given some notations in the section I would like to speak about the 
bases of scalar integrals in section [3] and to show that the classical base is not efficient, therefore I give 
a better base. Then, in section [H I would like to give a way to find coefficients in front of triangles in 
this new bases. Finally I finish by an exemple in section [5] the four-photon amplitudes. For a massive 
loop more work are needed. 



2 Notations 

So in this first section we gave some notations. The purpose of this article is to study the decom- 
position of a one-massless loop Feynman diagram amplitude. So consider a one-loop diagram with N 
ingoing legs (Fig. [I]), we note it amplitude: 



with Q the loop momentum in n = 4 — 2e dimensions and the denominator Df = Q\ + iX = 
(Q + Ti) 2 + iX. We decide to note the n-dimensional vectors in capital letters and the 4-dimensional 
vectors in small letters. The — 2e part of the loop momentum is fi. As the four and the — 2e parts 
Minkowski space are orthogonal, therefore the denominator "i" is written: 



"df" is the four dimensional part of the denominator "i". The function Num(Q) depends on the theory 
with described the loop. In classical gauge theory, this function is polynomial. 

In particularly we note I 7 ^ the amplitude of a one-loop diagram, with N external legs, where all 
internal propagators are scalar: 




(1) 



D f = Q? + iA = (qt + nf+iX = qf + iX-fi 2 = d? - /i 2 . 



(2) 




(3) 
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We often call those amplitudes: scalar integrals. There are well-known and we can express them with 
explicit analytic expressions. We recall some of them in Appendix [XJ Moreover, we use oftentimes the 
kinematical matrix <Sy and the Gram matrix G^, defined (for a massless loop) by: 

Sij = (qi - qj f Gij = 2pi-Pj (4) 




3 The bases of decomposition 



Now we assume a one-massless-loop Feynman diagram, described in a gauge theory by the amplitude 
(fTj) . The numerator of the amplitude is a polynomial function. This last point out is very importante 
to decompose the amplitude. 



3.1 The classical base 



As the loop are described in a standard gauge theory, therefore, the integrand of the amplitude is a 
rational function. We can expand it automatically in partial fractions. Each one gives a scalar integral. 
Therefore, we can write straightforward an amplitude like a linear combination of scalar integrals in n 
dimensions: 



N 



An = J]] calf + Rationnal terms + O(e). 



(5) 



Obviously the decomposition of an N external-leg one-loop Feynman diagram uses scalar integrals with 
at most N external legs. However, all those integrals {/", n E [1..N]} are not free. Indeed, we have the 
linear relation [TU [T5] : 



A' 



I N (S) = £ S^I^iS - - (-1)^+1 (N - n - 1) ^^^ +2 (S) 



det (S) 



(6) 



where S — {i} is the kinematical matrix obtained by eliminate the column and the line number "i" 
Moreover, we can show that the Gram determinant for N > 5 is equal to zero. So this last equation ([S 
simplify: 



N 



(7) 



det (S) 
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So we see easily that in the leading order in e, the reduction © simplify (as the loop is massless, we 
don't need tadpoles to decompose the loop): 

4 

A N = ^ ai I? + 0(e). (8) 

We need only scalar boxes, triangles, bubbles and rational terms to generate all the one-massless-loop 
Feynman diagrams. Generally we call the scalar function part: the analytic part and the rational terms 
the rational part. The call "analytic part" comes from fact that they contains poly logarithm functions 
of Mandelstam variables. The rational terms are studied in [TBIE]- But we are going to see, that this 
base B = {i^ 1 , .../J } is not efficient. 




The acknowledge of the coefficients {a^, ■■■ji} and the rational terms is enough to reconstruction all the 
amplitude. This base B = {Ig, ■■■12} 1S the most obvious, the canonical base, but not the more efficient. 
It is the subject of the next subsection. 



3.2 The problems of this base B 

The base of decomposition B = {1%, ---I^} nas a ^ least t w0 problems. 

The first comes from divergences. Indeed, the four-point functions are constituted by finite terms and 
infrared divergent terms, the three-point functions are infrared divergent and the bubbles: UV divergent. 
So three and four-point functions could have infrared divergences. Therefore, with the base B, we don't 
have separate explicitly the infrared, ultraviolet and finite structure. It could forget some compensations 
and it is not easy to found them. For instance, consider a loop without infrared divergences, therefore 
we can have three and four-point functions and not explicite compensations, this implies some numerical 
instabilities. 

The second problem concerns the famous Gram determinant. The reduction gives some negative 
powers of Gram determinants in front of the four-point functions. This determinant is spurious, because 
it could be equal to zero but this divergence is not physic, it has no physical explanation. Indeed, there 
is some unexplicit compensations to eliminate those determinants. 

Those two problems give the base B not efficient. So we have to find an other base. 
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3.3 A better base 



The main problem is the blend of infrared, ultraviolet and finite parts and in particularly, it should 
be genius if we could separate the infrared and the finite part of the four-point scalar functions. We 
could do it just by the application of the formula © to the four-point scalar integrals: 

II (S) = J2 ( S ~ «)< " (! + 2e ) Srl 7 ^ 2 ( 5 ) ' ( 9 ) 

The four-point scalar functions in n + 2 dimensions are totally finite and well-known, we recall them 
in Appendix [X] Moreover © see us that it appears one Gram determinant in the numerator of the 
coefficients in front of the four-point functions. So with (J5J), we reduce one time the problem of the Gram 
determinant, but we don't eliminate all the problem. So we conclude that the base B' = {if, i£,.Z£ +2 } 
is better than B. It remains to calculate the coefficients in front of each scalar integral. 




The ultraviolet divergences are carried out by the scalar bubbles, whereas the infrared structures are 
contained by the one and two-external-mass scalar triangles, but the finite structure are the four-point 
functions and the three-external-mass scalar triangle. So this decomposition partitions the amplitudes 
over the three kind of analytic structure without covering. 

We will see that this partition is very interesting. It eliminates many subtil compensations. And the 
result that we are going to show is that if the loop has no infrared divergences, so the decomposition 
has no infrared structure. The extension of this result to the ultraviolet one is not so easy. 



4 Coefficients in front of boxes and triangles in the base B = 

Now we are going to calculate the coefficients {a i , in this new base B' = {if, 7 3 n , /" + 2 }. In 

this section, we give the method to calculate the coefficients of the four and three-point functions. The 
ones in front of bubbles will be given in an other paper. 



4.1 Coefficients of four-point functions. 
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To transform the base B = {1% , I3 1 , 14 } to the base B = Ig, 1% +2 }, we have just to apply 

the formula © to the first base. So keeping only the four-point scalar function part, the transformation 
gives us: 



So, we obtain the bounding equation: 



where G (respectively S) are the Gram matrix (resp. kinematical matrix) of the four-point scalar 
function. The coefficient is given directly by the unitarity-cut method [7J [18]. We recall it, for 
example, consider a loop with N ingoing legs described by the amplitude (JTJ) and, we want the coefficient 
in front of the box obtained by pinching the denominators of the initial loop except the four ones called 
"a, b, c" and "<i". Therefore, the unitary cuts give us directly the coefficient by the limit : 

Num (g) d\d\dld\ 



N 



where df represents the four-dimensional part of the denominator D 2 : d 2 = qf + iX (eq. [2]). To obtain 
this limit we have just to solve the system {i £ [a..d], d 2 = 0}. The loop momentum is a four-dimensional 
vector, so we write it like a linear combination of four four-dimensional vectors (for example external 
ingoing legs). The system is linear, with four equations and four variables (four parameters of the linear 
combination of the loop momentum) , so we can solve it exactly, and we explicite the momentum of the 
loop. We note go this momentum. This process is explained in [5], I7|, [T2" l [T9 l 120] . 



4.2 Coefficients of three-point functions. 

In this subsection we use the Forde results and the Forde formalism given in [TT] . Here we give the 
method to find the coefficient in front the three-point functions, assuming that the four-point functions 
are in n + 2 dimensions. To simplify the proof, we use an amplitude with four external legs: 

Num (Q) 

which we decompose on the base B = {1%, I3 , /4" 1 " 2 } and rational terms: 

det (G) 
'det (S) 



M = I dnQ DlDlDiDV < 13 > 



Ai = -a ^f" ^ I2 +2 + J2 TlI h + J2 ^ J 2 + rationnal terms. (14) 



i=l i=l 



With the last subsection, we can compute, very easily the coefficient in front of the four-point function, 
we know a = Num(q ), where qo is the solve of the linear system {Vz £ [1..4], d 2 = given by the four 
cuts. Now we want to calculate the coefficients in front triangles. We are going to use the fact that 

det(G) +2 f fl( .x 
det (S) h ~ J d Q D 2 D 2 D 2 D 2 ' 

where &, = Ei=i ^ij 1 - The amplitude becomes: 



^4 = / ^"Q n2n2n2n2 -Num(q ) d n Q n2n2n2n2 



Num (Q) 



DiDiDiDi ™J - DiDiDiDj 

4 r d 2 2 

+ V T j; / d n Q ' + V jJZ + Rationnal Terms. (16) 

i=1 J L> l L> 2 L> Z U i f=1 
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For example, we assume that we want to calculate the coefficient in front of the three-point scalar 
function obtained by pinching the propagator number 1. So we apply, in the last equation (|16p . the 
linear application "Disc2,3,4", which cuts the three propagators 2,3 and 4 in four dimensions: 



Disc 2:3 ,4^4 = Disc 2 , 3: 4 J d n Q Jpg^^ = Num (g ) Disc 2 , 3 ,4 J d n Q i) , Z) , /);;/) , 

4 r n 2 

+ J2 T * Disc 2 , 3 ,4 / 

i=l ^ 1 



2()2n2n2 ' 



(17) 



We keep on only the three and four-point functions. The application "Disc2,3,4" inputs the three 
propagators d 2 ,, d 2 and d\ on-shell. As we keep only the four-dimensional part in the numerator, because 
we want the coefficient in front of the scalar integrals. Indeed, the — 2e dimensional part of the loop 
momentum in the numerator give the rational terms |16j . Therefore (|17j) becomes: 

rQ Num(Q) g ( ^ ^ 4) = ^ y d n Q ^-hD 2 g ^ ^ A)+Tl j rQ s ^ ^ 4) ( (18) 

where 6 (2, 3, 4) = 6 (d$) S (dj) S (df) . In this step we use the Forde Formalism [11]. As we have only 
three cuts, we have only three equations. We write again the loop momentum like a linear combination 
of four four-dimensional vectors {Kf,Kf, (K^K^), (K^ K$)} , where Kf and Kf are light-like 
vectors defined in the Appendix [B] But as we have only three equations, therefore, we can explicit only 
three parameters over four of the linear combination of the loop momentum. The one, which it remains, 
is noted "c". According to [11] . becomes: 



q 3 (c) = a Qi Kf + a 03 Kf + °- {K^K\) + ^± {K h A -fK\) . (19) 



where all a^- and vectors Ki are explained in Forde paper [11] and recalled in the Appendix[B] Moreover, 
we use the spinor notations introduced in [21] . Now we change the variables of integration and do the 
integrations over the three delta functions. We obtain: 

dc Jc = Num (q ) [ dc j J—^M +Tl [dc J c . (20) 

D\{c.y J D{{c) J 

If we want to calculate the coefficient of the three point function, we have just to solve the equation: 

l-Mf(co) - 0, (21) 
where D\ = d\ — /i 2 (eq. [2]). We put those / solves in the equation (|20|) and we obtain the result. 

Resultat 4.1 In the base B , the coefficient in front of the scalar integrals 1^ is: 

i 

Ti - h^Num^) , (22) 
»=i 

where cq is the solve of the equation ([UJ. For the other triangles, we have just to permute the cut 
propagators and the parameter b{. 



Here we see one of the interest of the base B rather than the base B. Indeed, not only the base B 
separate the infrared, ultraviolet and the finite parts but also the coefficients in front of scalar triangles 
are very simple to computed. There are no longer problems to obtain the free parameters, it is given by 
the equation (f2Tj) . The solve of this equation is often obvious. In the last example, if one of the legs p% 
or p4 is massless therefore, the solve is zero: Co = 0. 

To improve the computation in this base B , we are going to give some rules. With them, we are 
going to know directly, without computation, the null coefficients in front of triangles. 
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4.3 Infrared Divergences 



Proposition 4.2 Consider the decomposition of an one-massless-loop Feynman diagram on the base 
B = {I2 , ^3 1 12 +2 } ■ ^ e assume that this diagram has no infrared (soft or collinear) divergence. 
Therefore the coefficients in front of the scalar one or two-external-massive-leg triangles are zero. 

Proof : We can give an analytical proof, but with some arguments of reduction, the proof is obvious. 
Consider a one-loop-massless diagram free of IR divergences. The reduction by standard methods gives 
some sub-diagrams by pinching propagators. But those reduction cannot create infrared divergences. 
Indeed the pinched propagator plus the two- adjacent legs give a massive-external leg, which eliminates all 
infrared divergences around it. After the reduction, we obtain three-point sub-diagrams and the "finite" 
part of four-point scalar integrals. Since the three-point sub-diagrams are free of IR divergences, they 
cannot be expressed in term of one mass/two mass three-point scalar integrals and so the coefficients in 
front of them are zero. □ 

We can improve this result and consider a diagram with infrared divergences. 

Proposition 4.3 Consider the decomposition of an one-massless-loop Feynman diagram, with a soft 
divergence on the propagator number "i", on the base B = {/J, 1$, 12 +2 } ■ The coefficients in front of 
the one-external and two- external-massive-leg triangles are zero except the one which corresponds to the 
one-external-massive-leg triangle whom the external-massive-leg is opposite to the propagator number 





Proposition 4.4 Consider the decomposition of a one-massless-loop Feynman diagram, with a collinear 
divergence on the external leg "i", on the base {/J \ ^3 > 12 +2 } ■ The coefficients in front the one-external 
and two-external-massive-leg triangles are zero except the ones which correspond to the triangles, whom 
the leg number "i" doesn't belong to an external mass: 




E fli 7 < 




Proof : The proofs are the same like the first proposition 14.21 The external mass regular the infrared 
divergences. Therefore all three-point sub-diagrams have a null coefficient except which one preserves 
the divergence. □ 



Remark: In a loop, only photons, gluons or scalars create some divergent propagators. Indeed a 
fermion propagator has a numerator and this numerator compensate all infrared divergences. 



If we decide to decompose a one-massless-loop Feynman diagram on the baseB' = {/ 2 n , 7 3 ™, 7 4 n+2 }, 
using this last remark we deduce straightforward the three-point scalar integrals which have a null or 
non-null coefficients. 



5 An Exemple: the four-photon amplitudes 

We are going to calculate rapidly the helicity amplitudes of the four-photon amplitudes in QED 
for example, and in QED, scalar QED and supersymmetric QED ^ = 1. All the result are given in 
[HI 022 ESI HE W\ HE HE 123 EI] , but we are going to found the results. 
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The four-photon amplitudes 71 + 72 + 73 + 74 —> in QED are null at tree order. The first non 
null order is the one-loop order. Therefore, the amplitude is six one-massless-loop diagrams with four 
photons ingoing in a loop of fermions. So the decomposition of the amplitude on the base B is: 

Aa = ( a ' !2 +2 (&) + b ' J 3 (<?) + cI 2 +rationnal terms. (23) 

<r(l, 2,3,4) 

However with the last subsection, as the loop is a fermion loop, therefore, there are no infrared divergence 
so b = 0. The decomposition becomes: 

A4 = ^ (a T'l +2 (<t) + cl% (afj + rationnal terms. (24) 

2,3,4) 

a is obtained by the equation (lll| : 

where "a" is given by unitary-cut. If all photons have a positive or a negative helicity therefore, at least 
two adjacent photons have the same helicity, and as the loop is massless therefore the coefficient is null. 
We explain it in [20] - We have the same argument if three photons have the same helicity: 

a(+ + ++) = a(- + ++) = 0. (26) 

If we have two positive-helicity photons and two negative-helicity photons therefore the coefficient 
is non null only if the helicities ingoing alternately in the loop. The non null coefficient corresponds to 
the scalar integrals I% +2 (1 + , 2 _ , 3 + , 47). The argument in this integrals gives the order in which the 
photon ingoing to the loop. We alternate the negative and the positive helicity photons. With some 
calculation with find directly: 

4 (12) [34] e + u 2 

ft h) = — e - — -- — ; 27 

V ; [12] (34) s V 1 

To compute the rational terms and the coefficients of bubbles we need other methods. In [24l [25l [26j 
EH [281 E3 [301 E], we find: 

M+ + ++) = 8 ^ 2 ^fff + (28) 

^- +++ ^ 8 - 2 MSry + ° (e) ' (29) 

^(12) [34] f t 2 +u\ n+2 , 1onAS t-u 



M--++) = -8*^^-^|l + — ^—^(1324) + — 0(e). (30) 

The definition of each scalar integrals are given in Appendix El In the last helicity, the coefficients in 
front of bubbles are non null. They carried out ultraviolet divergences, but the amplitude is free of those 
divergence. But the soustraction of two bubbles compensate and those divergence disappear. We can 
extend those results to the six-photon amplitudes. All the results are given in [3171 131j . 



6 Conclusion 



In this paper, we suggest a way to decompose a one-massless-loop Feynman diagram. 

As the other methods of reduction, we decompose an amplitude like a linear combination of scalar 
integrals. But contrary to the methods already existing, here we use the scalar boxes in n+2 dimensions. 
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Thanks to this transformation, we win a degree of freedom which simplify the computation of the 
coefficients in front of the scalar triangles. 



But it remains works to finish the method. Indeed, we have to incorpore the computation of the 
coefficients in front of bubbles and rational terms. But as the example of the four-photon amplitudes 
shows, the bound between the amplitude and the ultraviolet divergences is not so simple as the one with 
the infrared divergences. 
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A Scalar integrals 

In this appendix, for sake of completeness, the definition of master massless integrals used in this paper 
is recalled, more details can be found in [32] • We also give det(G) the determinant of the Gram matrix 
Gij = 2pi.pj built with the external four momentum and det(5') the determinant of the kinematical 
S-matrix defined by iSy = (qj — qi) 2 where the qi are the four dimensional momentum flowing in the 
propagators. All was explain in section but we remember the draw. Moreover, we are going to use 



In [33], we find many formulae using those dilogaritms. We have to be careful, the Mandelstam variables 
could be negative. They are in the arguments of the dilogarithm, so to solve the problem of the Riemann 
sheet of the logarithm we use the analytic continuation by adding a small imaginary part, the prescription 
is: 




Figure 2: General structure of a loop. 



the dilogarithm function defined by [33] : 




(A.l) 



s — > s + iA, with A > 0. 



(A.2) 



And, all the scalar integrals must be multiplied by the angular integral: 



r r = 



r(i + e )r(i- e ) 2 

T(l - 2e) 



(A.3) 
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A.l Two-point Function. 



Pi x /^\/P4 

P2 -A_y^p 3 



S = S12 = S34 



5 2 = 



s 
s 



The determinants are given by: 



det(S 2 ) = s 2 
det(G 2 ) = s. 

The two-point function in n dimensions is: 

15 (s) = * (- S r = - e -\n(-s) + 2 + 0(e) 
and in n + 2 dimensions : 

/2 " +2(s) = ~2e(l-2e)(3-2 £ ) ( ~ s)1 



A. 2 One-external-massive-leg three-point Function. 

/* ~x Pi 

-^J S = S 3 4 

Pi 



/0 
5 3 ,i = 
\s 0/ 

The determinants are given by: 

det (5 3 ,i) = 
det (G 3 ,i) =-s 2 . 

This one-external-massive-leg three-point function in n dimensions is: 



and in n + 2 dimensions 



^ 2 < S » " „(! -.1(1 -a.) '-'' 
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A. 3 The two-external-massive-leg three-point Function 

P2 

Pi W- Pi f ml = 8X2 



Pi x V_y x I = s 34 







ml 


m 


^3,2 — j 












\m\ 









The determinants are given by: 

det(5 a>2 ) = 

det (G 3 , 2 ) = - (ml - mfj . 

The two-external-massive-leg three-point function in n dimensions is: 

r(m 2 m 2A _ 1 (-mjp-(-mj)- 
i 3 \m 1 ,m 2 ) — „ „ 



A. 4 The three-external-massive-leg three-point Function 



pi 



A 

P2 PS 




^3,3 = I ml 

1 ™3 

The determinants are given by: 

det(G3,3) = m\m\ — (mi.m 2 ) = - 
det(«S3 i3 ) = 2m\m\m\. 
The three-external-massive- leg three-point function in n dimensions is [32 : 



A 

¥ 



# (ml ml mf) = -L I f 2 Li 2 (l - — ) + 2Li 



ifWfiUW^ ■ 1 1 2 ' X1 





7T 2 




+ y 




-In 2 







where: 



rn 



Xl,2 = 



i+ml 



2 V V2/1 / V2/2/ V2/1 / \2/2 
± Va 



2m 2 



2/1,2 



2 — m\ + ml ± \/A 



2m 2 



A = + m 2 + w 3 — 1m\m\ — 2m\m\ — 2m\m\ — i sign(m 2 ) e. 
The formula (|A.21j) is available every where thanks to the small imaginary part 

v r K±ie, A > 



VA± 



le 



iiV^A, A < 
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A. 5 The zero-external-massive-leg four-point Function 





The determinants are given by: 

dot (G 4 , ) 
det (<S 4 ,o) 



S4.0 



/0 s 0\ 

t 



t 0/ 



s 



-2st(s + t) = 2stu 
t 2 s 2 = (24342) 2 . 



The zero-external-massive-leg four-point function in n dimensions is P_2_j : 

^o(s,t) = |i {(-*)-« + (-*)"«} -|f„(M), 

where: 

And in n + 2 dimensions, it is: 



A. 6 The one-external-massive-leg four-point Function 



/'1 



Pe 
/'1 



S = Sl2 
t = S 2 3 
U = S13 



S456 



/ 

s 



m\ s 0\ 

t 



t 0/ 



The determinants are given by: 

-2st (s + t- m 2 ) 
{stf = (lm3ml) 2 . 

The one-external-massive-leg four-point function in n dimensions is 



det(G 4 ,i) 
det(S 4 ,i) 



2stu 



1 



I2 A ( S ,t,m 2 ) = — {((-«)-« + (-*)-) + ((-«)"< - (-m 2 )-) + ((-*)" 



sie 2 
2 



Fi (s,t,m 2 ) 
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where: 



F 1 (s,i,m 2 ) =Li 2 ( 1- — ) +Li 2 ( 1 - — j -Li 2 y /y _1_ 



F (s,i) + <jLi 2 — 
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And in n + 2 dimensions, it is: 



7 4 "f (s,t,m 2 ) 



(-?) 


- Li 2 


to \ 




~J 


4} 


,to 2 ) 




-3) 





A. 7 The two-adjacent-external-massive-leg four-point Function 



P3 * K 
Pi 1 Pi 



t 



■ Su 
S425 



U = Si2 5 

m\ = s 25 
ml = s 3 6 



f 


2 

mf 


s 


°\ 


2 





2 


t 


s 


ml 










t 





0/ 



<?4,2A — 



The determinants are given by: 

det(G4 j2 yt) = —2s (m\m\ — t{m\ + to 2 , — s — £)) = — 2s(1toi4to 2 1) 
det(5 4 ,2A) = (si) 2 

The two-adjacent-external-massive-leg four-point Function in n dimensions is: 



!a,2A (s,t,ml,ml) 
where: 

F 2 a (s, t, to 2 , to 2 ,) = Li 2 ( 1 



1 



(sty 



{(- s )- £ + - (-to 2 )-) + ((-r e - (- 



— F 2 a (s,t,ml,ml) , 



Li 2 1 



5 ta (i) 



In 



In 



And in n + 2, it is: 



JTJa (s,t, mf ,m|) 



(fu — to 2 to!) (n — 3) 



F 2j4 (s,t, TO 2 ,TO 2 ) 



2to 2 to 2 + * (s - to 2 - TOl) 2 2 2 

; ^ — ox K (TOT, TO,, Too I 

2(n - 3) (tu - to 2 to 2 ) 3 v 2 ' 37 



A. 8 The two-opposite-external-massive-leg four-point Function 




S = Si4 3 
t = S 2 43 
W = S 2 3 
TO 2 = S14 

ml = s 56 
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( 


2 

m\ 


s 




2 

m\ 








t 


s 








m 2 . 




i 




0/ 



<$4,2B 



The determinants are given by: 

det(G4 i 2s) = — 2 (m\m\ — st) (m\ + mf — s — i) = 2it (si — m^m 2 ) 
det(S 4 , 2 s) = (st - to 2 ™ 2 .) 2 = (2toi3toi2) 2 = (2m 2 3m 2 2) 2 . 
The two-opposite-external-massive-leg four-point function in n dimensions is 

1 



7^25 (s,t,m 2 ,ml) = 



(st — m\m\)e 2 
1 

+ 



{((_„)-« _ (-m 2 )-) + ((-«)- - (-m 2 )-)} 
{((-<)- - (-to 2 )-) + ((-t)- - (-to 2 )-)} 



(si — m 2 m 2 )e 2 

— — 2 2 2 J2B (s, t, to 2 , m|) 
si — tojto^ 



where: 



Fib (s,t,m\, 



= - Li 2 1 - 



si 



Li 2 1 ! 



Li 2 1-^ + Li 2 1-21 + Li 2 1-^ +iln 2 f 



Fi (s,i,TO 2 ) +F X (s,t,m 2 ) - F (s,t) - S.U 2 (l - ~ y } • 



And in n + 2 dimensions, it is: 



rn+2 / 1 2 2 \ 



F 2 B(s,t 1 ml,ml 
u(n — 3) 



A. 9 The three-external-massive-leg four-point Function 



As 



Pi 



PA 



Ps 



/>(, 



Pi 



S = S1234 
t = S 7 i2 
U = S1256 
m| = S12 

m\ = s 34 
m§ = s 56 



^ 


to 2 


s 




m\ 





^2 




s 


m| 


Q 




V° 


t 







'4,3 



The determinant is: 

det(<S4,3) = (st — m\m^ 2 = (7mim2mi7) 2 = (7m 3 m2m 3 7) 2 . 
The three-external-massive-leg four-point function in n dimensions is: 

1 



11 (s,t,ml,ml,rnf) 



(si-m 2 m 2 )e 2 ^ ^ 



+ (-t)-'-(-ml)- e -{-rr4r 



2 

st — m\m\ 



F 3 (s,t,mf,ml,ml) , 
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where: 



F 3 (s,t,ml,ml,ml) 



= - 1]d | — r 1 In ( — ) - -In ( — I In I — 

+ Li 2 [ 1 1 - 



2 



Li 2 ( 1 - — ] - Li 2 ( 1 - 



A. 10 The four-external-massive-leg four-point Function 



P3 t P2 

Pi i p 7 

P6 



S = S1234 
t = S7812 



u = 
rni 



ml 
ml 



S1256 
= S12 
= S34 
= «56 
S78 



/ 




s 




mf 





ml 




s 


ml 







\ml 




m| 


0/ 



The determinant is given by: 

det(<S4.4) = (si — m\m\ — tttJto^) 2 — ^m\m\m\m\. 

The four-external-massive-leg four-point Function in n dimensions is: 

jn ( ± 2 2 2 2\ 
7 4 (s,t,m 1 ,m 2 ,m 3 ,m 4 ) 

where: 

17. / , 2 2 2 2\ 

F 4 (s,t,m 1 ,m 2 ,m 3 ,m 4 ) = 



and: 

P = ^/l-2A 1 -2A 2 + A?-2A 1 A 2 + A2 
m\m 3 

1 = 2 2 

(m 4 + ml) (ml + m\) 

? 9 
m 2 m 4 

2 = 2 2 * 

(mf + m^) (ml + m 3 ) 



— , ^2-^ (s,t,m 1 ,m 2 ,m 3 ,m 4 ) 

(mf + m 2 ) (m 2 +m^) p 



-Li 2 



Li 2 



1 - Ai + A 2 + /? 



1 - Ai - A 2 - p 
2T[ 



-Li 2 



Li 2 



1 — Ai + A 2 — p 



1 — Ai — A 2 + p 



2Ai 



1, f\l\, fl + \l-\2+P 

-2 ln ( v A|J ln ( v l + A 1 -A 2 -p 
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B Recall of the formalism introduced by Forde [llj. 



We have: 



and 



with 7 = p 3 .p4 ± \l {p3-Pif - p\p\- 



K; 



«03 



P3 - (p|/t) P4 

i - {pIpIH 2 ) 

Pi - (gI/7) P3 

i - (pIpIIi 2 ) ' 



pI (7 - Pi) 
7 2 - p\p\ 

p\ (7 - pf) 



(B.63) 



(B.64) 



r 



2 2 ' 

P3P2 
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